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Abstract We present a new model for the Galactic population of close double 
white dwarfs. The model accounts for the suggestion of the avoidance 
of a substantial spiral-in during mass transfer between a giant and a 
main-sequence star of comparable mass and for detailed cooling models. 
It agrees well with the observations of the local sample of white dwarfs 
if the initial binary fraction is ~ 50 % and an ad hoc assumption is made 
that white dwarfs with M ^ 0.3 Mq cool faster than the models suggest. 
About 1000 white dwarfs with V ^ 15 have to be surveyed for detection 
of a pair which has Mi + M2 Mch and will merge within 10 Gyr. 



1. INTRODUCTION 

The interest in the close double white dwarfs (hereafter CDWD) stems 
from several reasons: (i) white dwarfs are the endpoints of stellar evolu- 
tion; (ii) CDWD experienced at least two stages of mass exchange and 
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Table 1. Known close double white dwarfs and subdwarfs with WD companions 



AT 


Name 


p {d) 


Mx/Mq 




1 


WD 0135-052 


1.556 


0.47 


0.52 


2 


WD 0136+768 


1.407 


0.44 


0.34 


3 


WD 0957-666 


0.061 


0.37 


0.32 


4 


WD 1101+364 


0.145 


0.31 


0.36 


5 


WD 1204+450 


1.603 


0.51 


0.51 


6 


WD 1704+481A 


0.145 


0.39 


0.56 


7 


WD 1022+050 


1.157 


0.35 




8 


WD 1202+608 


1.493 


0.40 




9 


WD 1241-010 


3.347 


0.31 




10 


WD 1317+453 


4.872 


0.33 




11 


WD 1713+332 


1.123 


0.38 




12 


WD 1824+040 


6.266 


0.39 




13 


WD 2032+188 


5.084 


0.36 




14 


WD 2331+290 


0.167 


0.39 




15 


KPD 0422+5421 


0.09 


0.51 


0.53 


16 


KPD 1930+2752 


0.095 


0.5 


0.97 



Objects 1-14 are CDWD, objects 15 and 16 are sdB stars with suspected white dwarf 
companions. Mi is the mass of the brighter component and M2 is the mass of the companion. 
See the text for references. 



thus provide an important tool for testing the evolution of binaries; (iii) 
merging double CO white dwarfs are a model for SNela; (iv) CDWD 
may be the most important contributors to the gravitational wave noise 
at v ^ 0.01 Hz, possibly burying signals of other sources. 

We review the data on the currently known CDWD and present results 
of the modelling of the population of white dwarfs, which involves several 
new aspects. The most important are the treatment of the first stage of 
mass loss without significant spiral- in (Nelemans et al., 2000a), the use 
of detailed models for the cooling of white dwarfs and the consideration 
of different star formation histories for the Galactic disc. 

2. OBSERVED CLOSE DOUBLE WHITE 
DWARFS 

Table 1 lists the currently known CDWD for which the orbital pe- 
riod and the mass of at least one component are measured (Maxted and 
Marsh, 1999; Maxted et al., 2000a). The accuracy of the mass deter- 
minations is certainly not better than ±0.05 M Q . With the exception 
of WD 0135—052 and 1204+450, the brighter components of the pairs 
are likely to be He white dwarfs, since their mass is lower than 0.46 M Q , 
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the limiting mass for the helium ignition in a degenerate stellar core 
(Sweigart et al., 1990). In the range of M ~ (0.35 - 0.45) M white 
dwarfs could have CO cores and thick He envelopes (Iben and Tutukov, 
1985). However, the probability of the formation of these so called "hy- 
brid" WD is 4 - 5 times lower than for helium WD with the same mass 
(Tutukov and Yungelson, 1993; Nelemans et al., 2000b). The binary 
dwarfs WD 0957-666, 1101+364, 1704-481A, and 2331+290 are ex- 
pected to merge because of angular momentum loss via gravitational 
wave radiation. If both components are He dwarfs, the merger may 
result in the formation of a nondegenerate star or a supernova-scale ex- 
plosion (Nomoto and Sugimoto, 1977). In the case of CO companions, 
the formation of an RCrB star is expected (Webbink, 1984; Iben et al., 
1996). 

We included in Tab. 1 also data on two sdB stars with white dwarf 
companions (Orosz and Wade, 1999; Maxted et al., 2000b). Subdwarfs 
are supposed to be helium-burning stars. In these particular systems 
central helium burning will be completed before components merge. 
This makes these binaries candidates for future merging CO+CO white 
dwarfs. Remarkably, for KPD 1930+2752 the total mass of the system 
exceeds Mch, making it a possible SNIa candidate! 

About 10 more WD and sdB stars with suspected close WD compan- 
ions are known (Marsh, 2000; Maxted et al., 2001). However, for these 
systems the orbital periods or the masses of the components are not yet 
determined. 

3. FORMATION OF HELIUM WHITE 
DWARFS 

Three CDWD are of special interest - WD 0136+768, 0957-666, and 
1101+364. In these systems the masses of both components suggest 
that they are helium dwarfs and, thus, descend from degenerate cores of 
low-mass (sub)giants. If we designate main-sequence stars as MS, red 
(sub)giants as RG, white dwarfs as WD, the stage of mass exchange, 
either stable or unstable, as RLOF, and refer by subscripts 1 and 2 to 
the initial primary and secondary star, the evolutionary sequence which 
results in the formation of a double helium WD may be described as fol- 
lows: MSi + MS 2 -> RGi+MS 2 -» RLOF! -> WD!+MS 2 ^WD!+RG 2 
— > RLOF2 — > WD1+WD2. If the mass transfer is unstable, the change 
of component separation a is usually calculated by balancing the bind- 
ing energy of the envelope of the mass-losing star with the change of the 
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orbital energy (Webbink, 1984): 

Mi (Mi - M c ) _ 
An 



= a 



"M C M 2 Mi M 2 " 
2 af 2 a; 



(1) 



Here M c is the mass of the core of the mass-losing star, t\ is its radius, 
subscripts i and / refer to the initial and final values of the orbital sepa- 
ration, a is the efficiency of the deposition of orbital energy into the com- 
mon envelope and A is a parameter which depends on the density distri- 
bution in the stellar envelope; the usual assumption is A = 0.5 (de Kool 
et al., 1987). It has to be noticed, however, that for stars more massive 
than 3 M the values of A for red giants may be significantly lower (Dewi 
and Tauris, 2000). It would be worthwhile to investigate also lower mass 
stars, which form most of the observed CDWD. 

Giants with degenerate helium cores obey a unique core-mass - radius 
relation (Refsdal and Weigert, 1970). Neglecting a slight dependence on 
the total mass of the star, this dependence (for solar chemical composi- 
tion objects) may expressed as (Iben and Tutukov, 1985): 



where radius of the giant R and mass of the core M c are in solar units. At 
the instant when the star fills its Roche lobe, the radius given by Eq. (2) 
is equal to the radius of the Roche lobe. If the mass loss is unstable, the 
mass of the core of the mass-losing star doesn't grow during this stage 
and the mass of the new-born white dwarf is equal to M c . 

Applying Eq. (2) together with restrictions on the masses of the pro- 
genitors of helium white dwarfs (0.8 - 2.3 M ), Nelemans et al. (2000a) 
reconstructed the evolution of WD 0136+768, 0957-666, and 1101+364. 
Their conclusions may be summarized as follows: (i) in the first episode 
of mass loss, when the companion of the Roche lobe filling star was a 
main-sequence star of a comparable mass, no substantial spiral-in oc- 
curred, see Fig. 1; (ii) in the second mass loss episode, which resulted 
in the formation of the currently brighter white dwarf, the separation 
of the components was strongly reduced. The deposition of the energy 
into common envelope in this episode was highly efficient: in Eq. (1) 
the product aA ^ 3. A note has to be made in relation with the latter 
statement: since Eq. (1) provides nothing more than an order of magni- 
tude estimate, a formal solution which gives a > 1 indicates that energy 
deposited into common envelope has to be comparable to the orbital 
energy of the binary. 

Since the first mass loss episode is neither stable mass transfer nor a 
spiral-in in a common envelope and the physical picture of the process 
is absent, Nelemans et al. (2000a) suggested to use in the population 



R « 10 3 - 5 M c 4 , 



(2) 
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Figure 1. Periods after the first phase 
of mass loss as a function of the mass 
of the secondary component at this time. 
Solid, dashed, and dotted lines are for 
WD 0957-666, 1101+364, and 0136+768 
respectively. The top three lines are peri- 
ods needed to explain the mass of the last 
formed white dwarf. The middle three 
lines give the maximum period if the first 
white dwarf would be formed by conser- 
vative mass exchange. The lower three 
lines give the periods for the case when 
the formation of the dwarf was accompa- 
nied by a spiral-in (Eq. (1) with a\ = 2). 



Figure 2. Dependence of the relative 
variation of the separation of the compo- 
nents cif /en on the fraction of the mass 
lost by the star for the cases of spiral-in in 
a common envelope and AML regulated 
variation of a. Pairs of solid, dashed, and 
dash-dotted curves correspond to the ini- 
tial mass ratio of 1.0, 0.5, and 0.2, respec- 
tively. The upper curve of every pair is 
for the "AML formalism" [Eq. (3)]. 



synthesis studies a simple linear equation for the angular momentum 
balance: 

Ji-J { = l J i-^, ( 3 ) 

where J\ and Jf are the angular momenta of the pre- and post-mass- 
transfer binary, respectively, AM is the amount of mass lost by the 
binary and M is the total initial mass of the system. The parameter 7 is 
adjusted by fitting the orbital periods and masses of the three abovemen- 
tioned helium CDWD. It follows that 1.4 £7 £1.8. As Fig. 2 illustrates, 
when Eq. (3) is applied, the separation between the components changes 
much less drastically, compared to the case when the "standard" com- 
mon envelope formalism [Eq. (1)] is applied. For typical values of the 
fractional mass of the core Mj/Mj both formalisms give similar a//oj 
when q decreases to ~ 0.2. In the actual calculations, for a given Mf /Mi, 
the larger of two values of aj/ai was used. 
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4. COOLING OF WHITE DWARFS AND 
OBSERVATIONAL SELECTION 

The observed sample of CDWD is biased. Some objects were selected 
for study since their low mass already suggested binarity. White dwarfs 
must be sufficiently bright for the mass determination and the measure- 
ment of the radial velocities. This suggests to compare a magnitude 
limited model sample of dwarfs with the observations. Hitherto, follow- 
ing Iben and Tutukov (1985), it was assumed in all population synthesis 
studies that all WD can be observed for 10 8 yr, unless close pairs merge 
in a shorter time. The actual cooling curves were never applied. How- 
ever, recent studies (Driebe et al., 1998; Hansen, 1999; Benvenuto and 
Althaus, 1999; Sarna et al., 2000) show, that helium WD cool more 
slowly than carbon-oxygen ones. This is due to the higher heat content 
of the helium WD (Hansen, 1999) and residual nuclear burning in the 
relatively thick hydrogen envelopes (first noticed by Webbink (1975)). 
In our basic model, we take cooling curves for CO dwarfs from Blocker 
(1995) and for He ones from Driebe et al. (1998, hereafter DSBH), see 
Fig. 3. The initial models for these curves are obtained by mimicking the 
mass loss by stars and therefore may be considered as the most realistic. 
However, as we show below, cooling times given by them probably can 
not be taken at face value and need revision downward. Mass loss in 
stellar wind and thermal flashes, which extinguishes hydrogen burning 
(Iben and Tutukov, 1986; Sarna et al., 2000; Althaus et al., 2000), may 
provide the necessary mechanism. 




Figure 3. Cooling curves for 0.179, 
0.300, 0.414 M white dwarfs from 
Driebe et al. (1998) and for 0.6 and 
0.8 M Q ones from Blocker (1995), from 
right to left. Straight lines are the fits to 
the curves used in the simulations. 



' 4.0 6.0 8.0 10.0 12.0 

log t (yr) 

In addition to the sufficient brightness of the components, CDWD 
must have such orbital periods that radial velocities would be large 
enough to be measured, but small enough not to be smeared out during 
the integration. Following estimates of Maxted and Marsh (1999), we 



Close Double White Dwarfs 



7 



model this selection effect assuming that CDWD with 0.15 hr < P or {, < 
8.5 day will be detected with 100% probability and that above 8.5 day 
the detection probability decreases linearly to at P or b = 35 day. 

Our other basic assumptions may be listed as follows: we use Miller 
and Scalo (1979) IMF, flat distributions over initial mass ratios of the 
components and the logarithm of the orbital separation, and a thermal 
distribution of orbital eccentricities. 

5. AN OVERVIEW OF THE MODELS 

For the modelling of the population of double white dwarfs we used 
the numerical code SeBa (Portegies Zwart and Verbunt, 1996) with mod- 
ifications for low- and intermediate mass stars described by Nelemans et 
al. (2000b). Some of the important points are discussed below. 

Since most progenitors of white dwarfs are low-mass stars, the Galac- 
tic star formation history influences the current birthrate of CDWD and 
the properties of their population. This factor was not studied before. 
For the present study we take a time-dependent star formation rate in 
the Galactic disk as 

SFR(t) = 15 exp(-t/r) M yr" 1 , (4) 

where r = 7 Gyr. Equation (4) gives current SFR of 3.6 M yr _1 , which 
is compatible with the observational estimates (Rana, 1991; van den 
Hoek and de Jong, 1997). With this equation, the amount of matter 
that has been turned into stars over the lifetime of the Galactic disk (10 
Gyr in our model) is ~ 8 x 10 10 M . It is higher than the current mass of 
the disk, since a part of this matter is returned to the ISM by supernovae 
and stellar winds. We also compute several models with constant SFR, 
to allow comparison with previous work (see Tab. 2). 
The distribution of stars in the Galactic disk is taken as 

p(R, z) = po e~ R/H sech(z//i) 2 pc~ 3 , (5) 

where we use H = 2.5 kpc (Sackett, 1997) and h = 200 pc. The age and 
mass dependence of h is neglected. 

Table 2 gives an overview of our assumptions and model results and 
a comparison with some models of other authors. Model A is our basic 
model with an exponential star formation rate in the disk [Eq. (4)], 
initial fraction of binaries equal to 50% (i.e. with 2/3 of all stars in 
binaries) and cooling of white dwarfs according to DSBH, but modified 
as described below, in Sec. 6. The 50% fraction of binaries is suggested 
as a lower limit to their actual occurrence by the studies of normal main- 
sequence stars (Abt, 1983; Duquennoy and Mayor, 1991). This model, 
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Table 2. Summary of models 



Mod. 


SFR 


bin 


WD 


V 




SNIa 


I WD 


CWD 


Pwd,Q 


^PN,© 






% 


10 9 


10~ 2 


10~ 2 


10~ 3 


10~ 3 


10 8 


10 -3 


io- 12 


A 


Exp 


50 


9.2 


4.8 


2.2 


3.2 


4.6 


2.5 


19 


2.3 


A' 


Exp 


50 




3.0 


2.4 


2.2 


11.0 


1.0 






B 


Exp 


100 




8.1 


3.6 


5.4 


8.8 


4.1 






C 


Cns 


50 


4.0 


3.2 


1.6 


3.4 


3.6 


1.2 


8.5 


1.9 


D 


Cns 


100 




5.3 


2.8 


5.8 


6.1 


1.9 






HAN 


Cns 


100 




2.9 


2.8 


2.6 


23.0 


0.9 






ITY 


Cns 


100 


4.0 


8.7 


1.7 


1.9 


8.5 


3.5 


8.3 


1.5 


OBS 












4±2 






4-f- 20 


3 



The columns list the identifiers of the models, type of star formation rate assumed for the 
model (exponential or constant), initial fraction of binaries, total Galactic number of WD, 
rates of formation and merger of CDWD per 100 yr, rate of merger of CDWD with M1 + M2 > 
M(jh per 1000 yr (SNela), rate of formation of Interacting double WD (AM CVn type stars) 
per 1000 yr, total number of Close double WD in the Galaxy, local density of all WD per 
pc 3 , local rate of formation of planetary nebulae per pc 3 per yr. HAN and ITY denote (Ibcn 
ct al., 1997) and (Han, 1998) models, OBS - observational data. 

as well as all our models presented in the Table, were calculated with 
aX = 2 in Eq. (1) and 7 = 1.75 in Eq. (3). Model A' is similar to 
the model A, but assumes that the first stage of mass loss is a common 
envelope described by Eq. (1) instead of Eq. (3). Model B is similar to 
the model A, but has initially all stars in binaries, while model C has a 
constant star formation rate and 50% binaries. Model D has a constant 
star formation rate and 100% binaries. Models C and D were normalised 
to SFR of 4M yr -1 , to allow comparison with the models of Iben et 
al. (1997, ITY) and Han (1998). The former model was recalculated by 
LRY for the disk age of 10 Gyr, the age assumed in this study. Note, 
that ITY assume that unstable mass loss always results in the formation 
of a common envelope, and their formulation of the equation for the 
energy balance is different from Eq. (1). 

Briefly, the comparison of models shows the following. Models with an 
exponential SFR compared to the models with a constant SFR (mod. A 
vs. mod. C) have a higher number of old stars and a higher mass of the 
disk. This gives higher birthrates of CDWD and AM CVn stars. The 
rate of mergers giving SNela is similar, since it is determined mainly by 
the SFR in the past ~ 300 Myr (Tutukov and Yungelson, 1994). 

Model A with Eq. (3) for the first mass loss episode gives less mergers 
than model A' with Eq. (1), since it produces wider pairs. For the same 
reason the occurrence of SNela is lower in model A'. Formation rate of 
interacting systems (AM CVn's) is higher in model A', because, if two 
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common envelope phases occur, the second-formed WD is typically less 
massive than companion; if such a pair is brought into contact due to 
angular momentum loss via gravitational wave radiation, unequal masses 
of components favour stable mass transfer (Nelemans et al., 2001). 

Model D has an IMF rather similar to Han's model, but treats the 
first mass loss episode differently , does not have companion reinforced 
stellar wind, and has a higher qA value. This results in relatively less 
mergers in the first stage of mass loss and in higher birthrate and total 
number of CDWD and higher occurrence rate of SNe la. 

The Iben et al. model differs from all other models by applying a 
different equation for the evolution in the common envelope, which is, 
in practice, equivalent to the usage of much higher value of a\. This 
results in less frequent mergers in both stages of mass loss, hence, a 
higher total number of CDWD. The total number of the Galactic WD 
in this model is lower, compared to our models. This is due to the ITY 
assumption that close and wide binaries obey different distributions over 
the mass ratios of the components: flat for close systems and ocg -2 ' 5 for 
wide systems (see Tutukov k, Yungelson (1993) for details). 

Our models share with the Han and ITY models the same assumptions 
on the initial distributions of close binaries over mass ratios of compo- 
nents and their separations and have rather similar initial IMF for the 
primary components. The variations of the birthrates and numbers of 
objects in different classes (within factor ~ 3) arise mainly from the dif- 
ferences in the treatment of mass loss and transfer, in initial-final mass 
relations and other more minor details of population synthesis codes. 
Since the assumptions in the different studies are generally in agreement 
with results obtained from the modelling of stellar evolution, Tab. 2 il- 
lustrates the limits of the accuracy of predictions by the state of the art 
population synthesis studies for binary stars. 

Observational data which may help to constrain the models are rather 
scarce and uncertain, e.g., the estimates of local space density of white 
dwarfs p w d,Q differ by a factor 5 (Knox et al., 1999; Festin, 1998; Oswalt 
et al., 1995; Ruiz and Takamiya, 1995). Our basic model A, as well as 
model C, complies with these observational limits. Model A gives the 
annual birthrate of planetary nebulae in somewhat better agreement 
with the observations (Pottasch, 1996) than model C. 

Model A, as well as the rest of the models in Tab. 2, agrees reasonably 
with the occurrence of SNe la in galaxies similar to the Milky Way (E. 
Cappellaro, this volume). The difference in SNe la rates between the 
models has to be attributed mainly to the different treatment of mass 
loss and to the different initial-final mass relations for the components 
of binaries. 
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Figure 4- The model population of CDWD as a function of the orbital period and 
mass of the brighter dwarf of the pair. Top left: distribution of the currently born 
CDWD in model A. Top right: "observed sample" (Vn m = 15), with cooling according 
to DSBH and Blocker (1995). Bottom right: the same sample for the modified DSBH 
cooling. Bottom left: the total Galactic population of CDWD younger than 100 Myr. 
Dots with error bars: observed CDWD. 



6. MODELS VS. OBSERVATIONS 

1. Orbital periods and masses of close double white dwarfs. 

The parameters which are determined for all known CDWD are the 
orbital period and the mass of the brighter component of the pair. We 
plot in Fig. 4 the P Q rb — m distributions of the occurrence rate for the 
currently born CDWD and for the simulated magnitude limited samples 
(Mim = 15) for the models with different cooling prescriptions. 




og P (days) 



Figure 5. Cumulative distribution 
over the periods. The solid line is for 
the model with reduced cooling time 
for low-mass WD. The dashed line is 
for DSBH cooling without modifica- 
tions and the dash-dotted line is for 
a model with constant "observability" 
time of 100 Myr. Open squares give 
distribution of observed CDWD, filled 
circles give the distribution including 
the sdB+wd binaries. 
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In general, the white dwarf observed as the brighter member of the 
pair, is the one that was formed last, but occasionally, it may be the 
one that was formed first, due to the effect of differential cooling, see 
Fig. 3. The top right panel of Fig. 4 shows that if DSBH cooling curves 
are taken at face value, the observed sample contains predominantly 
low mass (M^O.3M ) white dwarfs, in contradiction to observations. 
However, as we mentioned above, low-mass WD may experience ther- 
mal flashes, which may reduce the mass of their hydrogen envelopes and 
extinguish hydrogen burning. This may be especially true for WD in 
close pairs (a^i?©), where one easily expects the formation of a com- 
mon envelope during a flash. Since estimates for the amount of mass 
which may be lost in a flash is not yet available , we make an ad hoc ex- 
treme assumption that white dwarfs with masses below 0.3 M Q cool like 
the most massive helium (0.46 M Q ) white dwarfs (hereafter we call this 
modified DSBH cooling). As can be seen from the bottom right panel 
of Fig. 4, this assumption brings the model in a much better agreement 
with observations. All model distributions which follow, are given for 
the modified DSBH cooling. 

For comparison, we plot in the left bottom panel of Fig. 4 the "ob- 
served" distribution assuming (like in the studies of other authors) that 
all WD are visible for 10 8 yr, unless a pair merges earlier due to GWR. 
Since in this case cooling curves are not applied, i.e. no magnitudes 
are computed for the white dwarfs, we can not construct a magnitude 
limited sample and this panel gives the total number of "potentially vis- 
ible" CDWD in the Galaxy. The better agreement with observations 
compared to the modified DSBH case is only apparent, as can be easily 
seen from the cumulative distributions (Fig. 5). 

Figure 5 shows a deficiency of observed systems between ~ 5 hr and 
~ 1 day. No selection effects are known that prevent detection of white 
dwarf binaries with such periods. This "gap" may be partially filled if 
we plot also subdwarf B stars with suspected white dwarf companions, 
thus assuming that the current sdB star is a white dwarf in the making. 
In addition to systems listed in Tab. 1, we include the binaries for which 
only P or b is determined. However, the number of detected systems is still 
too small to decide whether the "gap" is real and whether revisions of 
the stellar evolution models or CDWD formation scenarios are required. 

The merger of white dwarfs is one of the models for SNela (see, e.g., 
Livio 1999 for a review). We estimate that one merger candidate with 
Mi + Mi > 1.4 Mq may be found in a WD sample complete to V ~ 15 
which contains ~ 200 CDWD among a total of ~ 1000 WD. In view of 
the data in Tab. 2, this estimate is probably uncertain within a factor 
of at least 2 3. 
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Figure 6. Model sample of CDWD as a function of the orbital period and mass 
ratio of components for Vu m = 15 and the ratio of the luminosities of the components 
< 5. Left panel is for model A with the first phase of mass transfer treated by Eq. (3), 
right panel is the same for the model A' with standard Eq. (1) with aX = 2 in both 
mass transfer phases. Dots with error bars are observed systems. 

2. Period - mass ratio distribution. Our treatment of the first 
phase of mass transfer between a giant and a main-sequence star, which 
doesn't result in a significant spiral-in, leads to a concentration of the 
mass ratios of the model systems around q ~ 1. In practice, q in the 
observed systems can only be determined if the ratio of the luminosities 
of the components is ^ 5 (Moran et al., 2000). The P or b — Q distribution 
for the theoretical magnitude limited sample which obeys this criterion, 
is shown in Fig. 6. In model A' (with a common envelope in the first 
mass transfer) CDWD have predominantly g~ (1.75 — 2), which is not 
observed and there are hardly any systems with q~l. 

3. Mass spectrum of white dwarfs. The left panel of Fig. 7 shows 
the spectrum of white dwarf masses in a sample limited by V = 15 
and based on model A. It includes white dwarfs in close pairs which 
are brighter than their companions and genuine single objects, white 
dwarfs which are components of the initially wide pairs, merger products, 
white dwarfs which became single due to disruption of binaries by SNe 
explosions. In the same Figure we plot the cumulative distribution for 
DA white dwarfs masses, as estimated by Bergeron et al. (1992) and 
Bragaglia et al. (1995) 1 . 

The expected fraction of CDWD among all "observed" WD in our 
preferred model A is ~ 23%, slightly higher than the upper limit of the 
range suggested by Maxted and Marsh (1999) for DA white dwarfs: 1.7 
to 19% with 95% confidence. Since lowering the initial fraction of bina- 



1 These masses may have to be increased by about 0.05 Mq, if one uses models of white dwarfs 
with thick hydrogen envelopes for the estimates (Napiwotzki et al., 1999). 
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Figure 7. Left panel - mass spectrum of all white dwarfs in model A with modified 
DSBH cooling. Right panel - the same for a model in which all helium white dwarfs 
cool like a ~ 0.414 M dwarf and all CO white dwarfs cool like a 0.605 M dwarf. In 
both panels the theoretical cumulative distribution is shown as a solid black line and 
the cumulative distribution of observed systems as a grey line. 

ries below 50% would contradict the observations, the high percentage 
of CDWD in the model sample may mean that, even after our modifica- 
tion, DSBH results overestimate cooling times of the lowest mass white 
dwarfs. The right panel of Fig. 7 shows results of a simple numerical ex- 
periment in which we assign to all helium white dwarfs the cooling curve 
of a 0.414 Mq dwarf from DSBH and a cooling curve of a 0.605 M dwarf 
(Blocker, 1995) to all CO white dwarfs. This gives 17% for the fraction 
of CDWD. Even so, the contribution of the lowest mass WD to the total 
model sample still seems to be overestimated. However, there may exist 
yet unknown selection effects which prevent their detection. 

Yet another problem is the deficiency of (0.45 — 0.50) M white dwarfs 
in the model sample. In this part of the theoretical mass distribution 
only hybrid white dwarfs, descending via case B mass transfer from the 
stars in a narrow range of ~ (2.5 — 5) M do occur (Iben and Tutukov, 
1985). Their formation rate is relatively low. On the other hand, cooling 
of these objects which have almost pure helium envelopes with mass 
~ (0.01 — 0.1) M was not yet studied and it's possible that the simple 
interpolation between cooling curves for CO and He white dwarfs is not 
appropriate and gives misleading results. 

7. CONCLUSIONS 

1. The standard Algol- type evolution and spiral-in in common en- 
velopes cannot explain the wide orbits of the progenitors of close double 
WD after the first stage of mass loss. In low-mass binaries with similar 
masses of components, the loss of the envelope of the initial primary 
probably doesn't cause a strong spiral-in. In the absence of a physical 
picture of this process, we suggest to describe the change of separation 
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of the components by an angular momentum loss law [Eq. (3)]. An 
attempt to model the parameters of three observed CDWD with two 
helium components results in 1.4 7 1.8. 

2. Modelling of the observed sample of close binary WD shows that 
low-mass WD have to cool faster than is suggested by the recently 
published cooling curves for helium WD with thick hydrogen envelopes 
(Driebe et al., 1998; Sarna et al., 2000). We suggest that an additional 
mass loss by WDs in winds or thermal flashes may strongly diminish the 
masses of envelopes and reduce the cooling times. 

3. A reasonable agreement of the population synthesis results with 
observations in the estimates of the local space density of WD, and the 
Porb — m, Porb — Q, and N(P or b) distributions for CDWD may be achieved 
if the parameter 7 in Eq. (3) is ~ 1.7 and it's assumed that WD with 
M<O.3M cool like the most massive helium WD. Further reduction 
of cooling times for the helium WD may even improve agreement. 

4. The model with an initial binary fraction of 50% (2/3 of stars are 
in binaries) fits the observations better than the model in which initially 
100% of stars are in binaries. 

5. We estimate that the detection of at least one WD pair with 
Mi + M2 > Mch may require a survey of ~ 1000 white dwarfs. But 
the uncertainties in the input data for the population synthesis studies 
and badly known selection effects, most probably make this estimate 
accurate only within a factor of ~ 2 4- 3. 
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